Shock Wave Emissions of a Sonoluminescing Bubble 
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A single bubble in water is excited by a standing ultrasound wave. At high intensity the bubble 
starts to emit light. Together with the emitted light pulse, a shock wave is generated in the liquid 
at collapse time. The time-dependent velocity of the outward-travelling shock is measured with an 
imaging technique. The pressure in the shock and in the bubble is shown to have a lower limit of 
5500 bars. Visualization of the shock and the bubble at different phases of the acoustic cycle reveals 
previously unobserved dynamics during stable and unstable sonoluminescence. 

PACS numbers: 78.60.Mq, 43.25.-|-y 42.65. Re 



When intense ultrasound sound is applied to water, 
bubbles appear in the liquid. Among the properties they 
exhibit is sound radiation and emission of photons [0. 
In a controlled experiment, a single bubble alone may 
be driven stably in a standing ultrasound field. Here, 
intense light pulses of very short duration may be ob- 
served. Since this discovery |^ experimental work on 
the so called single bubble sonoluminescence (SBSL) has 
been extensively carried out to explain the phenomenon 
and the interesting features it displays: The energy is fo- 
cussed by 12 orders of magnitude [||, the light pulses are 
of picosecond duration |^ , the emitted light energy per 
pulse is in the MeV range, the blackbody-like spectrum 
peaks in the ultraviolet. The inter pulse synchronicity 
can be accurate on the picosecond scale ||] or chaotic 
on the microsecond scale Parameter studies have 
been done showing the region of stable SBSL lying on 
the boundary of a dissolution island |^ . Advanced driv- 
ing of the bubble is employed to increase the light output 
0. Theoretical and numerical work |8 -10| has been done 
to explain SBSL but so far few basic assumptions of the 
different theories could be verified experimentally. An 
inner shock wave launched in the interior of the bubble 
upon collapse has theoretically been assumed to account 
for the observed short SBSL light pulse and its spectrum 

Our experimental and numerical work focuses on the 
observation of shock waves being emitted into the liquid 
at bubble collapse |13|. The shock waves are visualized, 
the velocity of the front as it travels outwards is measured 
and the peak pressure of the shock is deduced. Effects ap- 
pearing at unstable SBSL are analyzed. The experiments 
are consistent with numerical simulations. In the experi- 
ment (figure 1^) the standing ultrasound wave is produced 
in a cylindrical cell filled with water of ambient temper- 
ature, distilled and de-gassed to 10-40% of ambient gas 
pressure. The cell consists of two piezoceramic cylinders 
connected by a glass tube [Q of 2.9 cm radius (overall 
height 12 cm). An optical glass plate closes the bottom, 
the top remains open. The driving frequency is 23.5 kHz 
and the driving amplitude « 1.2 to 1.5 bars. A bubble is 
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FIG. 1. Experimental setup 

inserted into the liquid with a syringe. The oscillating 
bubble is illuminated from the top by a copper vapour 
laser by light pulses of 7 ns duration (FWHV) followed by 
a low intensity tail of 30 ns. The wavelength is 511 nm. 
The repetition rate of one half the driving frequency is 
adjusted via a controllable delay to accommodate either 
locking to the driving signal or controlled phase shifting. 

Because shock waves modulate the phase of the laser 
light, optical filtering is used to transform this informa- 
tion into intensity modulations. Therefore the bubble im- 
age is passed through a (magnifying) 4/ spatial Fourier 
filter. Specifically, a Dark-Ground method is used 
that removes the zeroth order in the Fourier plane with 
a thin metal stick. Subsequently the image is picked up 
by a video camera delivering 25 frames/s. The shutter 
opening time is 0.25 ms such that the average image of 
2-3 shock waves is seen. Because of the stable repetitive 
bubble collapse a slow motion video of the oscillations 
]l^ , p^ and the shedding of shock waves is produced by 
slightly detuning the laser flash frequency. The images 
of the shock waves are digitized in a computer and their 
radius/time curves can be plotted. Because of their sub- 
micron size, sonoluminescing bubbles are hard to detect 
at collapse. But by recording the center of the shock wave 
the bubble position can be determined. Figure ^ shows 
the images of shock waves at different times. The shock 
is emitted at the main collapse of the bubble. The shock 
front shows up as a circle (figure ||a) and no anisotropy 
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FIG. 2. Images of shock waves emitted by a sonoluminesc- 
ing bubble, a) at t = 480 ns after collapse b) reflection from 
the side walls ai t = 3.54 fis before collapse. The arrow marks 
the bubble position; image side length is 3.5 mm. 

is seen within the optical resolution limit of 1.5 ^m, which 
is an indication of a symmetrical collapse. The front 
proceeds to the outer glass wall of the cylinder, reflects 
and moves inward again. The reflected shock wave has 
a duration > 40 ns and is distorted, presumably due to 
imperfections or misalignment of the glass wall. Figure 
shows the shock wave at the time it is refocused the 
most (3.54 /is before the next collapse). The main pres- 
sure peak seems to be « 700 /im away from the bubble 
at the lower end of a line structure. The refocused shock 
is sometimes powerful enough to kick the bubble through 
space a bit as it passes it. Weaker secondary reflections 
are also observable. At no time we could see a pressure 
pulse due to bubble rebounds The duration of the 
shock pulse can be determined to be 10 ns (FWHV). As 
this is on the order of the optical pulse length of 7 ns, 
this value is an upper bound. 

From successive images the velocity of the shock front 
is calculated. Figure ^ shows the average velocities as 
a function of distance from the bubble center. At very 
small distances (6-73 /im) an average value oi v — 2000 
m/s is measured, at larger distances the velocity of the 
shock front is decreasing to the ambient sound speed. 
Because the velocity of the shock decreases rapidly, the 
instantaneous velocity may be well above 2000 m/s. The 
pressure p in the shock can be determined from its ve- 
locity u by a Rankine-Hugeniot relation and a state 
equation for water, namely the Tait equation 
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FIG. 3. Average velocities {tave ~ 34 ns) of the 
SBSL-shock wave from successive images as a function of the 
distance from the generation (circles). The solid line is the 
mean velocity extrapolated by averaging over all measurement 
points up to a respective distance r from the bubble center. 
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p and p are the maximal density and pressure in the 
shock, po=998.2 kg/m'^ and po=l bar are the ambient 
density and pressure, 7i=7.025, 5=3046 bars Using 
(|l|), the shock pressure can be calculated to be 5500 bars. 

Numerical calculations have been carried out to further 
analyze the time dependence of the velocity and pressure 
of the shock front. The Gilmore model |2^ describing 
the radial motion of a bubble is used. 



RR 




C = 



\r=R 



p{R, R) 








p{R, R) + 




r=R 


= Co 1^ 


Po + B 














2a 











C C dt ' ^ ' 



(n-l)/2n 



4^ 
R 



R 



R is the bubble radius, C, p, and p are the speed of sound 
in the liquid, its density, and the pressure at the bub- 



ble wall, respectively. H — J 
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'^dp is the enthalpy 



of the liquid. Parameters were set to co=1483 m/s, 
(7=0.0725 N/m, /x=0.001 Ns/m^- a=i?o/8.86 is a hard- 
core van der Waals-term and k=5/3 the adiabatic 
exponent for argon [p^ . The pressure at infinity is 
Poo = Po + PeCos{2TT ft) , Pe and / are the driving pressure 
and frequency. 

The dynamics of the pressure pulse in the liquid is 
calculated by using the Kirkwood-Bethe hypothesis : 
the invariant quantity Y = R[H + ^R^) propagates with 
the characteristic velocity c + u, the local sound plus par- 
ticle velocity in the liquid. The outgoing characteristics 
are determined by l22] 
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Solving the bubble equation (||) gives the initial values i?, 
R, H and u = R ior each characteristic. Crossing char- 
acteristics in r — i space imply the generation of a shock. 
The exact position of the shock front can be obtained 
by equalization of the particle velocities in the hysteretic 
u — t curves |Q . 

Figure ^ shows the calculated shock wave velocity as 
a function of the distance from the bubble center. The 
maximal velocity of the pressure peak of 
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FIG. 4. Numerical calculations of a shock wave generated 
at collapse time of a bubble of 5 ^m ambient radius driven at 
23.5 kHz and 1.45 bars. Shown are data for the pressure peak 
that is travelling away from the bubble into the liquid as a 
function of distance from the bubble center. Dashed line: ve- 
locity of the peak, solid line: extrapolated mean shock veloc- 
ity, squares: particle velocity; the inset shows the calculated 
peak pressure in the shock as a function of distance from the 
bubble (solid line). The dotted line shows a r~ reference line 
for comparison. 

approx. 8300 m/s decreases within the first hundred /im 
to the ambient sound velocity. For comparison with the 
experiment the mean velocity of the peak is shown in 
figure ^. The experimentally obtained short time average 
and mean velocities compare quite well to the numerical 
findings. The particle velocity in the model reaches a 
maximum value of 333 m/s. The inset in figure ^ shows 
the peak pressure of the shock as it travels away from the 
center. The maximum value of 73000 bars at 1 /im decays 
quickly with increasing distance, within a hundred /im 
with a faster decay rate than the usual r~^. Though these 
numbers may be somewhat overestimated due to model 
limitations, it is seen that within the first few /im extreme 
conditions exist in the fluid. The greater dissipation close 
to the bubble may account for differences between our 
experimental results and previous inferences of the shock 
pressure from direct hydrophone measurements, which 
have yielded smaller values for the pressure . 

Using the shock wave as a microscope for the bubble 
position at collapse time, the time dependency and the 
position of the collapse have been measured for unsta- 
ble SBSL. Unstable SBSL occurs at the upper parameter 
values of the driving pressure and ambient gas concentra- 
tion: The ambient bubble radius grows until the bubble 
dynamics reaches an instability where bubble volume is 
rapidly lost. This cycle repeats itself on a slow time scale 
1^,^. Using rare events of double exposure at split-off 
time, the distance of the centers of two shock waves rep- 
resenting a bubble before and after the split-off can be 
used to calculate a lower bound of the bubble velocity 
due to the recoil of 0.5 m/s. During unstable SBSL the 




5 6 7 8 9 10 11 12 13 14 15 
time [s] 



FIG. 5. Measured radii of shock waves and position of the 
center during unstable SBSL. Data have been digitized from 
images taken each 40 ms at a constant phase of the driving, 
a: The upper line shows the radii of shock waves as a func- 
tion of time, as they change on a slow time scale. The lower 
two lines are relative x- and y-coordinates of the bubble at 
collapse, b: Zoom into the first seconds of upper figure. Spa- 
tial oscillations are seen as the shock radii (and the ambient 
radius of the bubble) are changing at e.g. t=8-9 s. 

bubble shows its dynamical behaviour over a long range 
of its ambient radius as a parameter. Figure |^ shows 
the radius of the shock wave and the bubble position at 
collapse time as a function of time. All experimental 
conditions were kept constant. It is seen that the bub- 
ble collapse does not occur at a constant phase any more. 
As the ambient bubble radius grows by diffusion, the col- 
lapse is shifted to later times. Because the illuminating 
flash occurs at a fixed phase of the driving signal some 
small time after the collapse, a later collapse decreases 
the time the shock front can travel outward until it is 
imaged. This way a larger ambient bubble radius shows 
up as a smaller shock radius. In figure |^a the recurrent 
process of growing on a slow time scale and a subsequent 
rapid decrease of ambient volume of the bubble is seen. 
At split-off the collapse time of the bubble is shifted by 
« 1 /ts with respect to the driving phase. Calculations of 
collapse time vs. bubble volume for SBSL-relevant bub- 
ble radii show (see also [^) that the bubble loses about 
one half of its volume. Most probably, multiple fragments 
(micro bubbles) will be generated. So far micro bubbles 



3 



have been observed at the lower amphtude threshold of 
SBSL, where they have a slowly decreasing velocity of 
initially 3-4 mm/s, do not return to the bigger bubble 
but dissolve within a few hundred fim within « 0.1 s. 
A closer look on a single cycle of figure |^ reveals, that 
the growing phase has a peculiar fine structure. Imme- 
diately after micro bubble split-off (e.g. at t=5.4 s, 10.4 
s) the collapse is shifted to later times (smaller shock 
radii), reaching a slightly decreasing plateau until it fi- 
nally increases. During each growing phase small bumps 
are seen (e.g. at ^=8-9 s). Looking at the position of 
the bubble one sees a connection: Each time the phase 
bumps, the bubble moves discretely through space and 
finally settles. The explanation may be oscillation in dif- 
ferent resonances, shock wave interaction or the acoustic 
field acting on the bubble is altered as it grows. 

We have visualized the generation of shock waves from 
a sonoluminescing bubble for the first time. Resulting 
from the enormous pressure inside the bubble and the 
great amount of energy transported by the surrounding 
liquid, the shock front is shown to have a faster speed 
than the ambient sound velocity. If one calculates p^ ] 
the change in the refractive index of water due to the 
theoretical local overpressure of 73 kbars, one arrives at 
a An = 0.23, which is 70% of the change at an air/water 
interface. Therefore measurements of the minimal radius 
by Mie scattering together with statements about the 
exact timing of the flash with respect to the minimal 
radius should be done keeping this in mind as the shock 
wave builds a scattering layer around the bubble. We 
cannot conclude from our data, whether the visualized 
shock in the liquid also consists of contributions of an 
hypothetical inner bubble shock that may be responsible 
for SBSL. The observed refocused shock can be shown 
to have an impact on the bubble dynamics. It would be 
interesting to see, if exact positioning and control of the 
timing of the reflected shock wave can be used to increase 
SBSL intensity @,§. 
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